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Based on concepts drawn from the ekpyrotic scenario and M-theory, we elab- 
orate our recent proposal of a cyclic model of the Universe. In this model, 
the Universe undergoes an endless sequence of cosmic epochs which begin 
with the Universe expanding from a 'big bang' and end with the Universe 
contracting to a 'big crunch.' Matching from 'big crunch' to 'big bang' is per- 
formed according to the prescription recently proposed with Khoury, Ovrut 
and Seiberg. The expansion part of the cycle includes a period of radiation 
and matter domination followed by an extended period of cosmic accelera- 
tion at low energies. The cosmic acceleration is crucial in establishing the 
flat and vacuous initial conditions required for ekpyrosis and for removing 
the entropy, black holes, and other debris produced in the preceding cycle. 
By restoring the Universe to the same vacuum state before each big crunch, 
the acceleration insures that the cycle can repeat and that the cyclic solution 
is an attractor. 

PACS number(s): 11.25.-w,04.50.+h, 98.80. Cq,98.80.-k 
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1 Introduction 



In a recent paper,§ we introduced the possibility of a cyclic Universe, a cos- 
mology in which the Universe undergoes a periodic sequence of expansion and 
contraction. Each cycle begins with a "big bang" and ends in a "big crunch," 
only to emerge in a big bang once again. The expansion phase of each cycle 
includes a period of radiation-, matter-, and quintessence-domination, the 
last phase of which corresponds to the current epoch of cosmic acceleration. 
The accelerated expansion phase dilutes by an exponential factor the en- 
tropy and the density of black holes and any other debris produced since the 
preceding big bang. The acceleration ultimately ends, and it is followed by 
a period of decelerating expansion and then contraction. At the transition 
from big crunch to big bang, matter and radiation are created, restoring the 
Universe to the high density required for a new big bang phase. 

Historically, cyclic models have been considered attractive because they 
avoid the issue of initial conditions.H Examples can be found in mythologies 
and philosophies dating back to the beginning of recorded history. Since 
the introduction of general relativity, though, various problems with the 
cyclic concept have emerged. In the 1930 's, Richard Tolmani discussed cyclic 
models consisting of a closed Universe with zero cosmological constant. He 
pointed out that entropy generated in one cycle would add to the entropy 
created in the next. Consequently, the maximal size of the Universe, and 
the duration of a cycle, increase from bounce to bounce. Extrapolating 
backwards, the duration of the bounce converges to zero in a finite time. 
Consequently, the problem of initial conditions remains. In the 1960's, the 
singularity theorems of Hawking and Penrose showed that a big crunch neces- 
sarily leads to a cosmic singularity where general relativity becomes invalid. 
Without a theory to replace general relativity in hand, considerations of 
whether time and space could exist before the big bang were discouraged. 
"Big bang" became synonymous with the beginning of space-time. However, 
there is nothing in the Hawking- Penrose singularity theorems to suggest that 
cyclic behavior is forbidden in an improved theory of gravity, such as string 
theory and M theory, and some people have continued to speculate on this 
possibility.!'! In the 1990's, observations showed that the matter density is 
significantly less than the critical density and that the scale factor of the 
Universe is accelerating.! Tolman's cyclic model based on a closed Universe 
is therefore observationally ruled out. 
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Curiously, the same observations that eliminate Tolman's cyclic model 
fit perfectly the novel kind of cyclic model proposed here. In our proposal, 
the Universe is flat, rather than closed. The transition from expansion to 
contraction is caused by introducing negative potential energy, rather than 
spatial curvature. Furthermore, the cyclic behavior depends in an essen- 
tial way on having a period of accelerated expansion after the radiation and 
matter-dominated phases. During the accelerated expansion phase, the Uni- 
verse approaches a nearly vacuous state, restoring very nearly identical local 
conditions as existed in the previous cycle prior to the contraction phase. 
Globally, the total entropy in the Universe grows from cycle to cycle, as Tol- 
man suggested. However, the entropy density, which is all any real observer 
would actually see, has perfect cyclic behavior with entropy density being 
created at each bounce, and subsequently being diluted to negligible levels 
before the next bounce. 

The linchpin of the cyclic picture is safe passage through the cosmic sin- 
gularity, the transition from the big crunch to big bang. In recent work with 
J. Khoury, B. Ovrut and N. Seiberg, we have proposed that a smooth tran- 
sition is possible in string theory.i'iln ordinary 4d general relativity, the big 
crunch is interpreted as the collapse and disappearance of four- dimensional 
space-time. Densities and curvatures diverge and and there is no sign that 
a transition is possible. But in the theory considered here, what appears to 
be a big crunch in the 4d effective theory actually corresponds to the mo- 
mentary collapse of an additional fifth dimension. As far as matter which 
couples to the higher dimensional metric is concerned, the three large spa- 
tial dimensions remain large and time continues smoothly. The temperature 
and density are finite as one approaches the crunch, and, furthermore, the 
geometry is flat just before and just after the bounce. In short, there is 
nothing to suggest that time comes to an end when the fifth spatial dimen- 
sion collapses. Quite the contrary, the most natural possibility is that time 
continues smoothly. Efforts are currently underway to establish this conclu- 
sion rigorously in string theory! The cyclic scenario considered here exploits 
this concept and is absolutely dependent on its validity. In the absence of 
a detailed theory of the transition from big crunch to big bang, we will pa- 
rameterize the bounce in terms of simple matching conditions incorporating 
energy and momentum conservation. 

The appeal of a cyclic model is that it provides a description of the history 
of the Universe which applies arbitrarily far back into our past. The model 
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presented here suggests novel answers to some of the most challenging issues 
in cosmology: How old is the Universe - finite or infinite? How large is 
it? What was the physical cause of its homogeneity, isotropy and flatness? 
What was the origin of the energy density inhomogeneities that seeded cosmic 
structure formation and are visible on the cosmic microwave sky? What is 
the resolution of the cosmic singularity puzzle? Was there time, and an 
arrow of time, before the big bang? In addition, our scenario has a number 
of surprising implications for other major puzzles such as the value of the 
cosmological constant, the relative densities of different forms of matter, and 
even for supersymmetry breaking. 

The cyclic model rests heavilv on ideas developed as part of the recently 
proposed "ekpyrotic Universe. "01 The basic physical notion is that the col- 
lision between two brane worlds approaching one another along an extra 
dimension would have literally generated a hot big bang. Although the orig- 
inal ekpvrosis paper focused on collisions between bulk branes and boundary 
branes,Q here the more relevant example is where the boundary branes col- 
lide, the extra dimension disappears momentarily and the branes then bounce 
apart.! The ekpyrotic scenario introduced several important concepts that 
serve as building blocks for the cyclic scenario: 

• boundary branes approaching one another (beginning from rest) corre- 
sponds to contraction in the effective 4d theoretic description;! 

• contraction produces a blue shift effect that converts gravitational en- 
ergy into brane kinetic energy;! 

• collision converts some fraction of brane kinetic energy into matter and 
radiation that can fuel the big bang;Bi! 

• the collision and bouncing apart of boundary branes corresponds to the 
transition from a big crunch to a big bangB 

A key element is added to obtain a cyclic Universe. The ekpyrotic scenario 
assumes that there is only one collision after which the interbrane potential 
becomes zero (perhaps due to changes in the gauge degrees of freedom on the 
branes that zero out the force). The cyclic model assumes instead that the 
interbrane potential is the same before and after collision. After the branes 
bounce and fly apart, the interbrane potential ultimately causes them to draw 
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together and collide again. To ensure cyclic behavior, we will show that the 
potential must vary from negative to positive values. (In the ekpyrotic ex- 
amples, the potentials are zero or negative for all inter-brane separations.) 
We propose that, at distances corresponding to the present-day separation 
between the branes, the inter-brane potential energy density should be posi- 
tive and correspond to the currently observed dark energy, providing roughly 
70% of the critical density today. That is, the dark energy that is causing 
the cosmic acceleration of the Universe today is, in this scenario, inter-brane 
potential energy. The dark energy and its associated cosmic acceleration play 
an essential role in restoring the Universe to a nearly vacuous state thereby 
allowing the cyclic solution to become an attractor. As the brane separation 
decreases, the interbrane potential becomes negative, as in the ekpyrotic sce- 
nario. As the branes approach one another, the scale factor of the Universe, 
in the conventional Einstein description, changes from expansion to contrac- 
tion. When the branes collide and bounce, matter and radiation are produced 
and there is a second reversal transforming contraction to expansion so a new 
cycle can begin. 

The central element in the cyclic scenario is a four dimensional scalar 
field 0, parameterizing the inter-brane distance or equivalently the size of 
the fifth dimension. The brane separation goes to zero as tends to — oo, 
and the maximum brane separation is attained at some finite value 4> ma x- 
For the most part our discussion will be framed completely within the four 
dimensional effective theory of gravity and matter coupled to the scalar field 
0. This description is universal in the sense that many higher dimensional 
brane models converge to the same four dimensional effective description in 
the limit of small brane separation. We shall not need to tie ourselves to 
a particular realization of the brane world idea, such as heterotic M theory 
for the purposes of this discussion, although such an underlying description 
is certainly required, both for actually deriving the scalar potential we shall 
simply postulate and for the ultimate quantum consistency of the theory. 
The extra dimensional, and string theoretic interpretation, is also crucial 
at the brane collision, where the effective four dimensional Einstein-frame 
description is singular and at which point we postulate a big crunch-big 
bang transition as outlined in Ref. (§). Again, for the present discussion we 
simply parameterize the outcome of this transition in terms of the density 
of radiation produced on the branes, and the change in the kinetic energy of 
the scalar field, corresponding to a change in the contraction/expansion rate 
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of the fifth dimension. 

The scalar field plays a crucial role in the cyclic scenario, in regularizing 
the Einstein-frame singularity. Matter and radiation on the brane couple to 
the Einstein frame scale factor a times a function /3(<f>) with exponential 
behavior as <p — > — oo, such that the product is generically finite at the brane 
collision, even though a = and <fi = — oo there. For finite 0, the coupling 
of the matter and radiation to <fi is more model-dependent. Models in which 
is massless at the current epoch, such as we describe in this paper, face 
a strong constraint due to the fact that <f) can mediate a 'fifth force', which 
is in general composition dependent and violates the equivalence principle. 
Again, without tying ourselves to a particular brane world scenario we shall 
consider models in which the coupling function /3(<f>) tends to a constant at 
current values of 4> (large brane separations), and the corresponding fifth 
force is weak. An example of such a model is the Randall-Sundrum model 
with the non-relativistic matter we are made of localized on the positive 



tension brane (see e.g. Ref. |10| for a recent discussion). In models where 
/3(0) does not tend to a constant at current values of (ft, one must invoke 
some physical mechanism to give the <fi field a small mass so that the fifth 
force is only short-ranged. This modification still allows for cyclic behavior, 
with an epoch of false vacuum domination followed by tunneling.0 

The outline of this paper is as follows. In Section ^|, we describe the req- 
uisite properties of the scalar field (inter-brane) potential and present a brief 
tour through one complete cosmic cycle. In subsequent sections, we focus in 
technical detail on various stages of the cycle: the bounce (Section [3]), pass- 
ing through the potential well after the big bang (Section f|), the radiation-, 
matter- and quintessence- dominated epochs (Section |), the onset of the 
contraction phase and the generation of density perturbations (Section ||). 
In Section [7|, we show that the cyclic solution is a stable attractor solution 
under classical and quantum fluctuations. In Section [8], we discuss the im- 
plications for the fundamental questions of cosmology introduced above. 



2 A Brief Tour of the Cyclic Universe 

The various stages of a cyclic model can be characterized in terms of a scalar 
field <p which moves back and forth in an effective potential V(<f)). In Sec- 
tion pTT| , we discuss the basic properties that V(</>) must have in order to 
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allow cyclic solutions. 

The stages of expansion and contraction can be described from two dif- 
ferent points of view. First, one can choose fields and coordinates so that 
the full extra-dimensional theory is reduced to an effective four- dimensional 
theory with a conventional Einstein action. The key parameters are the scale 
factor a and the modulus scalar field that determines the distance between 
branes. In this picture, the terms "big bang" and "big crunch" seem well- 
merited. The scale factor collapses to zero at the big crunch, bounces, and 
grows again after the big bang. However, what is novel is the presence of 
the scalar field <fi which is runs to — oo at the bounce with diverging kinetic 
energy. The scalar field acts as a fifth force, modifying in an essential way 
the behavior of matter and energy at the big crunch. Namely, the tempera- 
ture and matter density remain finite at the bounce because the usual blue 
shift effect during contraction is compensated by the fifth force effect due 
to 0. The arrangement seems rather magical if one is unaware that the 4d 
theory is derived from a higher dimensional picture in which this behavior 
has a clear geometrical interpretation. Nevertheless, for most of this paper 
we shall keep to the four dimensional Einstein description, switching to the 
higher dimensional picture only when necessary to understand the bounce, 
or to discuss global issues where matching one cycle to the next is important. 
The description of a single cycle from the 4d effective theory point-of-view is 



given in Section 2.2 



The same evolution appears to be quite different to observers on the vis- 
ible brane who detect matter and radiation confined to three spatial dimen- 
sions. In this picture, depending on the details, the brane is either always, or 
nearly always expanding except for tiny jags near the big crunch/big bang 
transition when it contracts by a modest amount. The branes stretch at 
a rate that depends on which form of energy dominates the energy density 
of the Universe. As the big crunch is approached, however, the expansion 
rate changes suddenly, and new matter and radiation is created (a brane has 
instantaneously collided into the visible brane and bounced from it). We 
describe some aspects of the visible brane viewpoint in Section |2.3| . 

This picture makes it clear that the big crunch does not correspond to 
the disappearance of all of space and the end of time but, rather, to the 
momentary disappearance of a fifth dimension. However, the behavior of 
gravity itself appears quite wild because it depends on the full bulk space- 
time, which is changing rapidly. One way of describing this picture is that 
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one has mapped the conventional big bang singularity onto the mildest form 
of singularity possible, namely the disappearance of a single dimension for 
an instant of time. Nevertheless there are delicate issues involved, as are dis- 
cussed in Ref. §, such as the fact that the effective four dimensional Planck 
mass hits zero at the singularity, so that gravitational fluctuations can be- 
come large. There are suggestions in specific calculations^ that physical 
quantities are nevertheless well behaved although a great deal more remains 
to be done to make the picture rigorous. 



2.1 The Effective Potential for a Cyclic Universe 

We will consider in this paper potentials V(<p) of the form shown in Figure 1, 
with the following key features: 

• The potential tends to zero rapidly as <fi — > — oo. One natural possibility 
for the extra dimension parameterized by is the eleventh dimension of M 
theory. In this case the string coupling constant g s oc e 7< ^, with some positive 
constant 7, and g s vanishes as —00. Non-perturbative potentials should 
vanish faster than any finite power of g s , i.e., faster than an exponential in 

• The potential is negative for intermediate 0, and rises with a region of 
large negative curvature, V"/V » 1 covering a range of of order unity 
in Planck mass units. This region is required for the production of scale 
invariant density perturbations, as proposed in Ref. [7| and detailed in Ref. |T2[ 
Attractive exponential potentials of this type could be produced, for example, 
by the virtual exchange of massive particles between the boundary branes. 

• As increases, the potential rises to a shallow plateau, with V"/V « 1 
and a positive height Vq given by the present vacuum energy of the Universe 
as inferred from cosmic acceleration and other astronomical evidence. The 
positive energy density is essential for having a cyclic solution since it pro- 
duces a period of cosmic acceleration that restores the Universe to a nearly 
vacuous state before the next bounce. The discussion here can be extended 
to potentials of a more general form. For example, it is not essential that 
the positive plateau persist to arbitrarily large since the cyclic solution 
only explores a finite range of > 0. Provided the condition V"/V <C 1 is 
satisfied over that range, the Universe undergoes cosmic acceleration when 
the field rolls down that portion of the potential. However, for simplicity, we 
will consider the example in Figure 1. 
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An explicit model for V(4>) which is convenient for analysis is 



\/(0) = K,(l-e- c V(0), (1) 

where, without loss of generality, we have shifted <fi so that the zero of the 
potential occurs at (f> = 0. The function F(<f>) is introduced to represent 
the vanishing of non-perturbative effects described above: F(<f>) turns off the 
potential rapidly as goes below 4> m i n , but it approaches one for > 4> m in- 
For example, F(<f>) might be proportional to e _1 / 9s or e" 1 ^ 93 , where g s oc e 7 ^ 
for 7 > 0. The constant Vo is set roughly equal to the vacuum energy 
observed in today's Universe, of order 1CT 120 in Planck units. We do not 
attempt to explain this number. Various suggestions as to how a suitable 
small positive vacuum energy could arise have been made.citesw,bp For large 
c, this potential has V"/V < 1 for <p k, 1 and V" /V > 1 for <p min < 
cf) < 0. These two regions account for cosmic acceleration and for ekpyrotic 
production of density perturbations, respectivelyBlia In the latter region, the 
constant term is irrelevant and V may be approximated by — Vo e~ C( ^ which 
may be studied using the scaling solution discussed in Section ^. 

For an arbitrary scalar potential of the form sketched, i.e. rising with 
negative curvature towards a flat plateau, the the scalar spectral index is 
given approximately by^ 2 "' 



n s ~ 1 — 4 



V \ 2 V"V 



V'J (V r 



(2) 



to be evaluated when the modes on the length scales of interest are generated, 
Stage (4) as described in Fig. 1. For the exponential form here, Eq. (Q) 
reduces to 

n s « 1 - -j. (3) 

Current observational limits from the cosmic microwave background and 
large scale structure data are safely satisfied for c = 10, which we shall 
adopt as our canonical value. 

The fact that the potential minimum is negative means that there are 
no strictly static solutions for (f> except anti-de Sitter space. However, as 
we shall show, the generic behavior - indeed an attractor - is a dynamical 
"hovering" solution in which roams back and forth in cyclic fashion between 
the plateau and — oo. The hovering solution is highly asymmetric in time. 
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Figure 1: The interbrane potential V((f>) versus 0, whose value (—00 < < 
0oo) determines the distance between branes. The shaded circle represents 
the maximum positive value of during the cycle. The various stages are: 
(1) quintessence/potential domination and cosmic acceleration (duration > 
trillion years); (2) kinetic energy becomes non-negligible, decelerated ex- 
pansion begins (duration ~ 1 billion years); (3) H = 0, contraction begins; 
(4) density fluctuations on observed scales created ((to^i?) 1 ^ 2 ~ 1 ms before 
big crunch); (5) kinetic energy domination begins {t m in ~ 1CT 30 s before 
big crunch); (6) bounce and reversal from big crunch to big bang; (7) end 
of kinetic energy domination, potential also contributes (t m i n ~ 10~ 30 s 
after big bang); (8) radiation dominated epoch begins tR ~ 1CT 25 s after big 
bang); (9) matter domination epoch begins (~ 10 10 s after big bang). As the 
potential begins to dominate and the Universe returns to stage (1), the field 
turns around and rolls back towards —00. 
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The field spends trillions of years or more on the plateau and mere instants 
traveling from the potential well to — oo and back. Gravity and the bounce 
provide transfers of gravitational to kinetic to matter-radiation density that 
keep the Universe forever hovering around the anti-de Sitter minimum rather 
than being trapped in it. 



To set the context for the later sections, we present a brief tour through a 
single cycle, using the labels in Figure 1 as the mileposts. Stage 1 represents 
the present epoch. The current value of the Hubble parameter is H — (15 
billion yr)^ 1 . We are presently at the time when the scalar field is acting as 
a form of quintessence in which its potential energy has begun to dominate 
over matter and radiation. Depending on the specific details of the potential, 
the field may have already reached its maximal value (grey circle), turned 
back, and begun to evolve towards negative values. If not, it will do so in 
the near future. Because the slope of the potential is very small, <fi rolls very 
slowly in the negative direction. As long as the potential energy dominates, 
the Universe undergoes exceedingly slow cosmic acceleration (compared to 
inflationary expansion), roughly doubling in size every Ho 1 = 15 billion 
years. If the acceleration lasts trillions of years or more (an easy constraint 
to satisfy), the entropy and black hole densities become negligibly small and 
the Universe is nearly vacuous. The Einstein equations become: 



where H is the Hubble parameter and G is Newton's constant. We will gen- 
erally choose 87rG = 1 except where otherwise noted. Accelerated expansion 
stops as V((f)) approaches zero and the scalar field kinetic energy becomes 
comparable to the potential energy, Stage (2). The Universe continues to 
expand and the kinetic energy of scalar field continues to red shift as the 
potential drops below zero. A nearly scale invariant spectrum of fluctuations 
on large length scales (beyond our current Hubble horizon) begins to develop 
as the field rolls down the exponentially decreasing part of the potential. 
The evolution and perturbation equations are the same as in the ekpyrotic 



2.2 The View from Effective 4d Theory 




(4) 




(5) 
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modelJ3li2l Solving these equations, one finds that the decelerated expansion 
continues for a time Hq 1 /c, which is about one billion years (c is the pa- 
rameter in V((f>), Eq. ([I])). At Stage (3), the potential becomes sufficiently 
negative that the total scalar field energy density hits zero. According to 
Eq. (fD, H = and the Universe is momentarily static. From Eq. (^), a < 0, 
so that a begins to contract. The Universe continues to satisfy the ekpyrotic 
conditions for creating density perturbations. Stage (4), about one second 
before the big crunch, is the regime where fluctuations on the current Hubble 
horizon scale are generated. As the field continues to roll towards — oo, the 
scale factor a contracts and the kinetic energy of the scalar field grows. That 
is, gravitational energy is converted to scalar field (brane) kinetic energy dur- 
ing this part of the cycle. Hence, the field races past the minimum of the 
potential at Stage (5) and off to — oo, with kinetic energy becoming increas- 
ingly dominant as the bounce approaches. The scalar field kinetic energy 
diverges as a tends to zero. At the bounce, Stage (6), matter and radiation 
are generated, the scalar field gets a kick and increases speed as it reverses 
direction, and the Universe is expanding. Through Stage (7), the scalar ki- 
netic energy density (oc 1/a 6 ) dominates over the radiation (oc 1/a 4 ) and the 
motion is almost exactly the time-reverse of the contraction phase between 
Stage (5) and the big crunch. As the field rolls uphill, however, the small 
kick given the scalar field and, subsequently, the radiation become impor- 
tant, breaking the time-reversal symmetry. The Universe becomes radiation 
dominated at Stage (8), at say 10~ 25 s after the big bang. The motion of <fi is 
rapidly damped so that it converges towards its maximal value and then very 
slowly creeps downhill. The damping continues during the matter dominated 
phase, which begins thousands of years later. The Universe undergoes the 
standard big bang evolution for the next 15 billion years, growing structure 
from the perturbations created when the scalar field was rolling downhill at 
Stage (4) . Then, the scalar field potential energy begins to dominate and cos- 
mic acceleration begins. Eventually, the scalar field rolls back across = 0. 
The energy density falls to zero and cosmic contraction begins. The scalar 
field rolls down the hill, density perturbations are generated and <fi runs off 
to — oo for the next bounce. The evolution in terms of conventional variables 



is summarized in Figure [L2 
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Figure 2: Schematic plot of the scale factor a(t), the modulus <p(t), and 

one cycle, where t is Einstein frame proper 

time. The scale factor starts out zero but expands as t^, and the scalar 
field grows logarithmically with t, in the scalar kinetic energy dominated 
early regime. Then, when radiation begins to dominate we have a oc t?, and 
the scalar field motion is strongly damped. This is followed by the matter 

2 

era, where a oc ts, and a potential dominated phase in which a(t) increases 
exponentially, before a final collapse on a timescale Hq 1 , to a = once more. 
H 5 is proportional to the proper (five dimensional) speed of contraction of 
the fifth dimension. To obtain a cyclic solution, the magnitude of at the 
start of the big bang, H 5 (out), must be slightly larger than the value at the 
end of the big crunch, H 5 (in). This is the case if more radiation is generated 
on the negative tension brane (see Appendix). 

13 



2.3 The View from the Visible Brane 



Thus far, we have described the evolution in terms of the usual Einstein 
frame variables, a and <fi. However, as emphasized in the next Section, these 
variables are singular at the transition from big bang to big crunch, and they 
do not present an accurate picture of what an observer composed of matter 
confined to the brane would actually see. As a approaches zero, the density of 
matter and radiation scale as l/(a/9(0)) 3 and l/(a(3(<f))) 4 , respectively, where 
(3 (<f>) is a function of <fi which scales as 1/a as a tends to zero. Therefore the 
densities of matter and radiation on the branes are actually finite at a = 0. 

This scaling of the density with af3(4>) rather than a can be understood 
rather simply. First, the spatial volume element on the branes is that induced 
from five dimensions. When the brane separation is small, one can use the 
usual formula for Kaluza-Klein theory, 

d4 = e-yft+ds\ + e 2 yft+dy 2 , (6) 

where ds\ is the four dimensional line element, y is the fifth spatial coordinate 
which runs from zero to L, and L is a parameter with the dimensions of 
length. If we write the four dimensional line element in conformal time 
coordinates, as ds 2 = a 2 (—dr 2 +dx 2 ), then since from the Friedmann equation 
we have (a' /a) 2 = |(0') 2 > we see that a is proportional to in the big 

crunch. Hence a three dimensional comoving volume element d 3 xa 3 e 
remains finite as a tends to zero. Thus the density of massive particles tends 
to a constant. What about the density of radiation? First, recall the usual 
argument that the energy of a photon diverges at a = 0. Consider a set of 
comoving massive particles in a spacetime with metric a 2 r]^ where r]^ is the 
Minkowski metric. The four velocities of the particles obey u^u u g^ v — — 1. 
Hence, if they are comoving (u = 0), then we must have u° = a" 1 . Now 
a photon moving in such a spacetime has a constant four- momentum, p M = 
E(l,n), with n 2 = 1. The energy of the photon, as seen by the comoving 
particles, is —u t *p u g lll/ = E/a, which diverges a tends to zero. However, 
in the present context, the metric to which the comoving particles couple 

is e~^/^^a 2 rj^ v . Therefore, we have u° = a~ 1 e^^ and the energy of the 
detected photons is finite as a tends to zero. In other words, the effect of the 
scalar field approaching — oo is precisely such as to cancel the gravitational 
blueshift. 
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The second crucial use of the higher dimensional metric is in piecing to- 
gether the global view of the spacetime. If one only had the Einstein frame 
scale factor a, it would not be clear how to match to the next cycle, since 
a = at the bounce. But the scale factor on a brane, a/3(<f>), is nonzero at 
each bounce and may be so matched. In fact, in the examples studied in this 
paper, the scale factors ao and a>i (which are the brane scale factors in the 
simplest models) both undergo a net exponential expansion within a cycle, 
and decrease for very brief periods - either just before the brane collision (for 
ao) or just after it (for a\). An observer on either brane would view the cos- 
mology as one of almost uninterrupted expansion, with successive episodes of 
radiation-, matter-, and quintessence-domination ending in a sudden release 
of matter and radiation. 

Both matter and radiation are suddenly created by the impact of the 
other brane. The forewarning of this catastrophic event would be that as 
(3 (<f>) started to rapidly change, one would see stronger and stronger violations 
of the equivalence principle (a 'fifth force'), and the masses and couplings of 
all particles would change. In the case of M theory, the running of the string 
coupling to zero would presumably destroy all bound states such as nucleons 
and send all particle masses to zero. 

3 Through the Bounce 

To have repeating cycles, the Universe must be able to pass smoothly from 
a big crunch to a big bang. Conventionally, the curvature and density sin- 
gularity when the scale factor a approaches zero has been regarded as an 
impassable obstacle to the understanding of what came 'before' the big bang. 
However, the brane world setup sheds new light on this problem. The key 
feature is that the apparent singularity in the effective four-dimensional de- 
scription corresponds to a higher dimensional setup in which the four dimen- 
sional metric is completely non-singular. When the extra dimension (or outer 
brane separation) shrinks to zero, there is no associated curvature singularity, 
and the density of matter on the branes remains finite. The most conserva- 
tive assumption, based on the higher dimensional picture, is that the branes 
bounce from (or, equivalently, pass through) each other and time continues 
smoothly, with some conversion of brane kinetic energy to entropy. The sep- 
aration of the two branes after the bounce corresponds to re-expansion in 
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the four- dimensional effective theory. 

How can this be reconciled with the singular four-dimensional descrip- 
tion? The point explained in Ref. |8] is that the usual four-dimensional vari- 
ables, the scale factor a and the scalar field 0, are a singular choice at a = 0. 
Each is poorly behaved as the branes collide, but in the brane picture physi- 
cal quantities depend on combinations of the two variables that remain well- 
behaved. These nonsingular variables may be treated as fundamental, and 
matching rules derived to parameterize the physics of inelastic brane colli- 
sions. If the system conjectured in Ref. || be properly embedded 
within string theory, the matching conditions will be derivable from funda- 
mental physics. 



3.1 Non-singular variables 

The action for a scalar field coupled to gravity and a set of fluids pi in a 
homogeneous, flat Universe, with line element ds 2 = a 2 {r)(—N 2 dT 2 + dx 2 ) is 



S = J d 3 xdr N- 1 (-3a' 2 + ^V 2 ) - N ((a/3) 4 Y, lPi + a*V (</>)) 



(7) 



We use t to represent conformal time and primes to represent derivatives 
with respect to r. N is the lapse function. The background solution for the 
scalar field is denoted 4>(t), and V(4>) is the scalar potential. 

The only unusual term in ([?[) is the coupling of the fluids pi, which we 
treat as perfect fluids coupled only through gravity. The action for a perfect 
fluid coupled to gravity is just — / d 4 x^—gp, where the density p is regarded 
as a function of the coordinates of the fluid particles and the spacetime met- 
ric.0 For a homogeneous isotropic fluid, the equation of state P(p) defines 
the functional dependence of p on the scale factor a, via energy-momentum 
conservation, dlnp/dlna = —3(1 + w), with w = P/p. For example, for 
radiation, p oc a -4 and for matter p oc a~ 3 . 

We assume these fluids live on one of the branes, so that rather than 
coupling to the Einstein-frame scale factor a, the particles they are composed 
of couple to a conformally related scale factor a/3(0), being the scale factor 
on the appropriate brane. For simplicity we have only written the action for 
fluids on one of the branes, the action for fluids on the other brane being a 
xerox copy but with the appropriate (3{4>). 
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The function /3(0) may generally be different for the two branes, and for 
different brane world setups. But as mentioned above there is an important 
universality at small separations corresponding to large negative 0. In this 
limit, which is relevant to the bounce, the bulk warp factor becomes irrelevant 
and one obtains j3 ~ e~^^, the standard Kaluza-Klein result. This behavior 
ensures that a(3 is finite at collision and so the matter and radiation densities 
are, as well. 

The equations of motion for gravity, the matter and scalar field are 
straightforwardly derived by varying (0) with respect to a, N and 0, after 
which N may be set equal to unity. Expressed in terms of proper time t, The 
Einstein equations are 



H 



2 8vrG/l- 2 



+ V + (3 4 p R + (3 /k p M ^ , (8) 



\ = ~ (> -V + /3W^Vm) , (9) 

where a dot is a proper time derivative. As an example, we consider the 
case where there is radiation and matter (pm) on the visible brane only, 
which could in principle be either the positive or negative tension brane. 
Then the above equations are supplemented by the dynamical equation for 
the evolution of 0, 

+ 3#0 = -V %4> - P^Pm (10) 
and the continuity equation, 

where a = a/3(0) and p is the pressure of the fluid component with energy 
density p. Note that only the matter density contributes to the 0-equation, 
because, if pr oc l/(a/?) 4 , the radiation term is actually just a constant 
times N in the action, contributing to the Friedmann constraint but not the 
dynamical equations of motion. 

If /3(0) is sufficiently flat near the current value of 0, these couplings 
have modest effects in the late Universe, and the successes of the standard 
cosmology are recovered. For example the total variation in since nucle- 
osynthesis is very modest. In Planck units, this is of order (tr/im) 1 where 
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t r is the time at the beginning of the radiation dominated epoch and nu- 
cleosynthesis begins at ijv ~ 1 sec. It is utterly negligible for values of t r 
earlier than the electroweak era. However, the coupling of matter to <fi pro- 
duces other potentially measurable effects including a 'fifth force' causing 
violations of the equivalence principle. Current constraints can be satisfied 
if M Pl (\n(3)^< 10' 3 EE 

As the Universe contracts towards the big crunch, a — > 0, the scalar field 
runs to — oo and the scalar potential becomes negligible. The Universe be- 
comes dominated by the scalar field kinetic energy density since it scales as 
a -6 whereas matter and radiation densities scale as a -3 and a -4 respectively 
(ignoring the (3 factor). As scalar kinetic domination occurs, the scale factor 
a begins to scale as (— £) 3 , and the background scalar field diverges logarith- 
mically in time. The energy density and Ricci scalar diverge as (— 1)~ 2 , so 
that t = is a 'big crunch' singularity. 

As explained in Ref. in the simplest treatment of brane world models 
there is only one scalar field modulus, the 'radion,' which runs off to minus 
infinity as the scale factor a approaches zero. The singular variables, a and 
<ft, can be replaced by the non-singular variables: 

a = 2acosh((0 — 0oo)/v6) ai = — 2asinh((0 — (p^/VQ). (12) 

The kinetic terms in the action define the metric on moduli space. In terms 
of the old variables one has the line element — 3da 2 + ^a 2 d(p 2 , and a = is 
clearly a singular point in these coordinates. However, in the new coordinates 
in Eq. (|12"1), the line element is |(— da 2 , + da 2 ), which is perfectly regular even 

when the Einstein frame scale factor a = \\J ai — a 2 vanishes, on the 'light- 
cone' an = ai. For branes in AdS, an and ai are the scale factors on the 
positive and negative tension branesQ so that (3 = 2cosh((0 — 0oo)/v6 or 
— 2sinh((0 — 0oo)/v^6) respectively for matter coupling to these branes. 

Notice that the constant field shift (p^ is arbitrary. Its effect is a Lorentz 
boost on the (a ,ai) moduli space. In the Kaluza-Klein picture (|6]), a con- 
stant shift in (f) can be removed by rescaling four dimensional spacetime 
coordinates and re-defining the length scale L of the extra dimensions. In 
the absence of matter which couples to 0, or of a potential V(0), this shift is 
unobservable, a reflection of the global symmetry — > 0+ constant of the 4d 
effective theory. However, this symmetry is broken by V(4>), and by matter 
couplings. In fact, the scale factor a\ must be positive in order for it to be 



18 



interpretable as a 'brane scale factor', and this requires that <p < 0oo- 

We shall find it convenient to choose = to be the zero of the potential 
V ((/)), and then to choose (p^ so that never vanishes for the solutions we 
are interested in. (In fact, since a\ has a positive kinetic term,§ a suitable 
coupling to moduli fields will always guarantee that a\ 'bounces' away from 
zero.ll! In this paper, for simplicity we ignore this complication by picking 
large enough that no such 'bounce' is necessary.) 

Both cio and a\ are 'scale factors' since they transform like a under rescal- 
ing space-time coordinates. However, unlike a they tend to finite constants 
as a tends to zero, implying an alternative metrical description which is not 
singular at the 'big crunch'. In the brane world models considered in Ref. ^|, 
cio and a% actually represent the scale factors of the positive and negative ten- 
sion branes respectively. Since there are no low energy configurations with 
cio < ai, the 'light cone' do = a\ is actually a boundary of moduli space and 
one requires a matching rule to determine what the trajectory of the system 
does at that point. A natural matching rule is to suppose that at low energies 
and in the absence of potentials or matter, the branes simply pass through 
one another (or, equivalently, bounce) with the intervening bulk briefly dis- 
appearing and then reappearing after collision. This rule was detailed in 
Ref. D, where simple models satisfying the string theory background equa- 
tions to all orders in a' were given. In the Appendix we discuss the collision 
between boundary branes in terms of energy and momentum conservation, 
and the Israel matching conditions. 

Let us now comment on the character of the trajectory in the (do, de- 
plane. The Friedmann constraint reads 

«'o 2 - a!l = ~ ((d/3) 4 p + i(og - a?) V(0 O )) • (13) 

If the energy density on the right hand side is positive, the trajectory is 
time-like. If the right hand side is zero (for example if the potential vanishes 
as 0o - > — oo and if there is no matter or radiation), then the trajectory is 
light-like. If the right hand side is negative, the trajectory is space-like. 
The trajectory for the cyclic solution in the do — dx plane is shown in 



Figure |3.1| . The insert shows a blow-up of the behavior at the bounce in which 
the trajectory is light-like at contraction to the big crunch (the Universe is 
empty) and time-like on expansion from the big bang (radiation is produced 
at the bounce). In these coordinates, the scale factor increases exponentially 
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over each cycle, but the next cycle is simply a rescaled version of the cycle 
before. A local observer measures physical quantities such as the Hubble 
constant or the deceleration parameter, which entail ratios of the scale factor 
and its derivatives in which the normalization of the scale factor cancels out. 
Hence, to local observers, each cycle appears to be identical to the one before. 

3.2 From Big Crunch to Big Bang 

In this section we solve the equations of motion immediately before and after 
the bounce, and discuss how the incoming and outgoing states are connected. 
The nonsingular 'brane' scale factors ao and a\ provide the natural setting 
for this discussion, since neither vanishes at the bounce. As emphasized 
above, the Einstein frame scale factor a, and the scalar field are singular 
coordinates on field space at the bounce. Nevertheless, since our intuition is 
much better in the Einstein frame, we shall also give formulae for a and <j) 
near the bounce. In subsequent sections, we shall frame the discussion almost 
entirely in terms of Einstein frame variables, for the most part using the 
nonsingular variables ao and a\ solely as a 'bridge' connecting the incoming 
big crunch to the outgoing big bang. 

Before the bounce there is little radiation present since it has been expo- 
nentially diluted in the preceding quintessence-dominated accelerating phase. 
Furthermore, the potential V(<f>) becomes negligible as runs off to minus 
infinity. The Friedmann constraint reads (a' /a) 2 = |0' 2 , and the scalar field 
equation, (a 2 0')' = 0> where primes denote conformal time derivatives. The 
general solution is 



that r < before collision. A is an integration constant which could be set 
to unity by rescaling space-time coordinates but it is convenient not to do 
so. The Hubble constants as defined in terms of the brane scale factors are 
a' /al and a^/aj which at r = take the values +\~ 3 H 5 (in) and —\~ 3 H 5 (in) 
respectively. 




(14) 



where A 



e <^ o/v / 6_ We choose r = to be the time when a vanishes so 
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Figure 3: Schematic plot of the a -ai plane showing a sequence of cycles 
of expansion and contraction (indicated by tick marks). The dashed line 
represents the "light-cone" a = a\ corresponding to a bounce (a = 0). Each 
cycle includes a moduli kinetic energy, radiation, matter and quintessence 
dominated phase and lasts an exponentially large number of e-folds. The 
insert shows the trajectory near the big crunch and bounce. The potential 
energy V{<p) assumed takes the form shown in Fig. 1. 
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Re-expressing the scalar field as a function of proper time t = J adr, we 
obtain 

<j> = y|ln (j# 5 (m)t) . (15) 

The integration constant H^{in) < has a natural physical interpretation as 
a measure of the contraction rate of the extra dimension (See equation (^)): 




where L 5 = LeV^ is the proper length of the extra dimension, L is a 
parameter with dimensions of length, and £5 is the proper time in the five- 
dimensional metric, 

dh = ae'^dr = e'^dt, (17) 

with t being FRW proper time. Notice that a shift 0oo can always be com- 
pensated for by a rescaling of L. As the extra dimension shrinks to zero, H§ 
tends to a constant, H 5 (in). 

If the extra dimension shrinks adiabatically and backreaction from parti- 
cle production can be ignored, then the matching rule conjectured in Ref. |8] 
states that H 5 after the bounce should be given by H 5 (out) = —H 5 (in). 
However, if radiation is produced, H 5 (out) takes a different value. If one 
is given the densities of radiation produced on both branes, then H 5 (out) 
may be inferred from energy and momentum conservation, and the Israel 
matching conditions, as we show in the Appendix. 

Immediately after the bounce, scalar kinetic energy dominates and 
remains nearly constant, as shown in Figure |2.2| . The kinetic energy of the 
scalar field scales as a -6 and radiation scales as a -4 , so the former dominates 
at small a. It is convenient to re-scale a so that it is unity at scalar kinetic 
energy-radiation equality, t r , and denote the corresponding Hubble constant 
H r . The Friedmann constraint in Eq. (|T3|) then reads 



2 

and the solution is 



a') 2 = l -H 2 r {l + a- 2 ), (18) 




3^ [2*TH£{out)H{ 

{H r T + 23) 




(19) 
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The brane scale factors are 



ao = a 



( A -V/V6 + Xe -4>/^) = A (x{l + ^f) + X-^Hhhou^r 



= A(x{l + ^f- )-2*\- 1 H?H£((wt)Ty (20) 

! 1 

Here the constant A = 2& (H r / H$(out)) 3 has been defined so that we match 
ao and a% to the incoming solution given in flT3]). As for the incoming solution, 
we can compute the Hubble constants on the two branes after collision. They 

5 - i 

are ±X~ 3 Hs(out) + 2~s X~ x Hr H§ on the positive and negative tension branes 
respectively. 

For H r < 2%X- 3 H 5 , the case of relatively little radiation production, 
immediately after collision ao is expanding but a x is contracting. Whereas 
for H r > 22X~ 3 Ha, both brane scale factors expand after collision. We 
shall concentrate on the former case in this paper, in which we are near the 
adiabatic limit. If no scalar potential V(4>) were present, the scalar field 
would continue to obey the solution (0), converging to 




2i H5M)\ (21) 



r 



This value is actually larger than 0^ for H r < H 5 X~ 3 2% , the case of weak 
production of radiation. However, the presence of the potential V(<p) alters 
the expression ( plj) for the final resting value of the scalar field. As crosses 
the potential well traveling in the positive direction, H 5 is reduced to a 
renormalized value H 5 (out) < H 5 (out), so that the final resting value of the 
scalar field can be smaller than 0oo. If this is the case, then a\ never crosses 
zero, instead reversing to expansion shortly after radiation dominance. (In 
the calculations of Ref [T3, where we assumed the potential vanished after 



collision, this effect did not occur. Instead, we invoked a coupling of a\ to a 
modulus field which caused it to bounce off a x = 0.). 

If radiation dominance occurs well after <ft has crossed the potential well, 
Eq. (|2l|) provides a reasonable estimate for the final resting value, if we use 
the corrected value H^(out). The dependence of ( pl| ) is simply understood: 
while the Universe is kinetic energy dominated, a grows as tz and <p increases 
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logarithmically with time. However, when the Universe becomes radiation- 
dominated and a oc £2, Hubble damping increases and converges to the 
finite limit above. 

4 Across the Well 

Using the potential described in Section ^7i] and, specifically, the example 
in Eq. (|l|), this Section considers the motion of back and forth across the 
potential well. We will show that evolution converges to a stable attractor 
solution. Our main purpose, though, is to explore the asymmetry in the 
behavior before and after the bounce that is an essential component of the 
cyclic solution. 

Over most of this region, V may be accurately approximated by — Vn e~ c< ^. 
For this pure exponential potential, there is a simple scaling solution^ 

a(t) = \t\^ V = -V e-^ = - P -^^, p=|, (22) 

which is an expanding or contracting Universe solution according to whether 
t is positive or negative. (We choose t = to be the bounce.) From the 
expression for V, we see that varies logarithmically with \t\. 

At the end of the expanding phase of the cyclic scenario, there is a period 
of accelerated expansion which makes the Universe empty, homogeneous and 
flat, followed by rolling down the potential V(<j>) into the well. After has 
rolled sufficiently and the scale factor has begun to contract (past Stage (3) in 
Figure |2.1| ), the Universe accurately follows the above scaling solution down 
the well until encounters the potential minimum (Stage (5) in Figure |2.1|) . 

Let us consider the behavior of under small shifts in the contracting 
phase. In the background scalar field equation and the Friedmann equation, 
we set = (pB + <50 and H = H B + 5H, where 0# and H B are the background 
quantities given from fl2"2"|). To linear order in 5(p, one obtains 

50 + ^^0-^50 = 0, (23) 

with two linearly independent solutions, 50 ~ and t 1_3p , where p< 1. In 
the contracting phase, the former solution grows as t tends to zero. However, 
this solution is simply an infinitesimal shift in the time to the big crunch: 
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5<p oc 0. Such a shift provides a solution to the Einstein-scalar equations 
because they are time translation invariant, but it is physically irrelevant 
since it can be removed by a redefinition of time. The second solution is a 
physical perturbation mode and it decays as t tends to zero. Hence, we find 
that the background solution is an attractor in the contracting phase. 

We next consider the incoming and outgoing collision velocity, which we 
have parameterized as H 5 (in) and H 5 (out) in the previous Section. Within 
the scaling solution (p2|), we can calculate the value of incoming velocity by 
treating the prefactor of the potential F(4>) in Eq. ([[]) as a Heaviside function 
which is unity for > min and zero for < 4> m in, where 4> m in is the value of 
at the minimum of the potential. We compute the velocity of the field as 
it approaches m j n and use energy conservation at the jump in V to infer the 
velocity after m j n is crossed. In the scaling solution, the total energy as 
approaches min from the right is |0 2 + V = 3p 2 /t 2 , and this must equal the 
total energy |0 2 evaluated for just to the left of m i n . Hence, we find that 
= \/6p/t = jQpVmin/ (1 — 3p) at the minimum and, according to Eq. (PL 
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At the bounce, this solution is matched to an expanding solution with 

H 5 (out) = -(1 + x)H s (in) > 0, (25) 

where x is a small parameter which arises because of the inelasticity of the 
collision. 

In order to obtain cyclic behavior, we shall need x t° be positive or, 
equivalently, the outgoing velocity to exceed the incoming velocity. There 
are at least two effects that can cause x t° be positive. First, as we discuss 
in the Appendix, x is generically positive if more radiation is generated on 
the negative tension brane than on the positive tension brane at collision. 
Secondly, x can g e t a positive contribution from the coupling of (3 ((f)) to the 
matter created on the branes by the collision; see Eq. (0). Both effects are 
equally good for our purposes. For the present discussion, we shall simply 
assume a small positive x is given, and follow the evolution forwards in time. 

Since x is small, the outgoing solution is very nearly the time reverse 
of the incoming solution as starts back across the potential well after the 
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bounce: the scaling solution given in fl2"2"|) , but with t positive. As time 
proceeds however, the contribution of x becomes increasingly significant. In 
the time-reversed scaling solution, H 5 tends to zero. For x > H 5 remains 
positive and overshoots the potential well. Vq is exponentially smaller than 
the kinetic energy density at the bounce, so even a tiny fraction x suffices to 
reach the plateau after crossing the potential well. 

We can analyze this overshoot by treating x as a perturbation and using 
the solution in Eq. (^3|) discussed above, <50 ~ t" 1 and t l ~ Zp . The latter is 
a decaying mode in the contracting phase before the bounce but it grows in 
the expanding phase. One can straightforwardly compute the perturbation 
in 5Hz, in this growing mode by matching at <p m in as before. One finds SH^ = 
YlxH® j(? where H§ is the background value, at the minimum. Beyond this 

point, 5H 5 grows as t^l c oc eV^, for large c, whereas in the background 
scaling solution H§ decays with as e^v^s - ^ 2 ^. When the perturbation is of 
order the background value, the trajectory departs from the scaling solution 
and the potential becomes irrelevant. The departure occurs when the scalar 
field has attained the value 

2 c 2 fl2x\ 2 

<l>Dep = 0mm+ -l n \ V \^\~^r) \ V min\- (26) 

As passes beyond (pDep the kinetic energy overwhelms the negative poten- 
tial and the field passes onto the plateau Vq with H 5 nearly constant (see 
Figure Ol) , and equal to 




Mout) «xhr H s(in), (27) 



until the radiation, matter and vacuum energy become significant and H§ is 
then damped away to zero. 

Before moving on to discuss these late stages, it is instructive to compare 
how rapidly travels over its range before and after the bounce. The time 
spent to the left of the potential well (0 < m j„) is essentially identical in 
the incoming and outgoing stages for x << 1, namely 



3\/2| V m in\ 



(28) 
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For the outgoing solution, when has left the scaling solution but before 
radiation domination, the definition Eq. (|T6|) may be integrated to give the 
time since the big bang at each value of </>, 



The time in Eq. (|29]) is a microphysical scale. The corresponding formula for 
the time before the big crunch is very different. In the scaling solution ( p2]) 
one has for large c 



e c(fli-(/) m i, l )/2 Q e c(if>-(f> min )/2 
K<P) = -\lTTTl = ~2 Kinl (30) 



The large exponential factor makes the time to the big crunch far longer than 
the time from the big bang, for each value of <fi. This effect is due to the 
increase in H 5 after the bounce, which, in turn, is due to the positive value 
of X- 



5 The Radiation, Matter and Quintessence 
Epochs 

As the scalar field passes beyond the potential well, it runs onto the positive 
plateau Vo- As mentioned in the last section, the value of H 5 (out) is nearly 
canceled in the passage across the potential well, and is reduced to H 5 given 
in Eq. (|27j). Once radiation domination begins, the field quickly converges 
to the large t (Hubble-damped) limit of Eq. (|19|) , namely 

c = y|ln(2i J H 5 (out)/ J ff r )), (31) 

where H r is the Hubble radius at kinetic-radiation equality. The dependence 
is obvious: the asymptotic value of depends on the ratio of H^(out) to 
H r . Increasing H^(out) pushes <fi further, likewise lowering H r delays radia- 
tion domination allowing the logarithmic growth of in the kinetic energy 
dominated phase to continue for longer. 
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As the kinetic energy red shifts away, the gently sloping potential gradu- 
ally becomes important, in acting to slow and ultimately reverse 0's motion. 
The solution of the scalar field equation is, after expanding Eq. flT9| ) for large 
r, converting to proper time t = J a{r)dT and matching, 



V3H, 



r - a~ 3 /* dtcfV^, (32) 



where as above we define a(t) to be unity at kinetic-radiation equal density. 
During the radiation and matter eras, the first term scales as t~i and t~ 2 re- 
spectively. For a slowly varying field, V<p is nearly constant, and the potential 
gradient term in Eq. (53) scales linearly with t, so it eventually dominates. 

When does <fi turn around? We give a rough discussion here, ignoring 
factors of order unity. First, we use the fact that V ~ t$ 2 where t is the 
present age of the Universe, and roughly we have to ~ 10 5 t m , where t m is 
the time of matter domination. As we shall see, <fi may reach its maximal 
value 4> m ax and turn around during the radiation, matter or quintessence 
dominated epoch. All three possibilities are acceptable phenomenologically, 
although the case where turnaround occurs in the radiation epoch appears 
more likely. For example, <p m ax is reached in the radiation era, if, from 
Eq. ©, " 

i ^ 

1 1() ; (— rf^OrjY 1. (33) 



If the Universe becomes radiation dominated at the GUT scale, t r ~ 10~ 25 
seconds. Then, only if we fine-tune such that V/V > 10 10 does t max exceed 
t m . This corresponds to the case where we have 10 10 e-foldings or more of 
cosmic acceleration at late times as <fi rolls back, far more than required for 
the cyclic solution. The bound changes somewhat if the Universe becomes 
radiation dominated as late as nucleosynthesis (t r ~ 1 second). In that 
case, even if V /V^(((pc) is not much greater than unity, the scalar field turns 
around in the matter era or later. For turnaround in the matter era, we 
require 

1 i 

3xl0_4 ^ (r)"(?; (0c) ) 3 ^ 30 - (34) 

Finally, if the field runs to very large 0c*, so that V^/V^c) ~ ce~ C(/>c is expo- 
nentially small, then <ft only turns around in the quintessence-dominated era. 
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For the example considered here, the natural range of parameters corresponds 
to turnaround occurring during the radiation-dominated epoch. Hence, by 
the present epoch, the field is rolling monotonically in the negative direction 
and slowly gaining in speed. Consequently, the ratio of the pressure to the 
energy density is increasing from its value at turnaround, w — — 1, towards 
zero. Depending on the details of the scalar potential V(0), it is conceivable 
that the increasing value of w could ultimately be observationally detectable. 

Once the field has turned around and started to roll back towards the 
potential well, the second term in Eq. (|32"D dominates. For our scenario 
to be viable, we require there to be a substantial epoch of vacuum energy 
domination (inflation) before the next big crunch. The number of e-foldings 
N e of inflation is given by the usual slow-roll formula, 



for our model potential. For example, if we demand that the number of 
baryons per Hubble radius be diluted to below unity before the next con- 
traction, which is certainly over-kill in guaranteeing that the cyclic solution 
is an attractor, we set e 3Ne ^ 10 80 , or N e ^ 60. This is easily fulfilled if 4>c 
is of order unity in Planck units. 

From the formulae given above we can also calculate the maximal value 
4>c in the cyclic solution: for large c and for t r » x~ 1 tmin, it is 



where we used H r 1 ~ t r , the beginning of the radiation-dominated epoch. 
From Eq. ( |36l) we obtain 



This equation provides a lower bound on t r . The extreme case is to take 
\V min \ ~ 1. Then using V ~ KT 120 , c ~ 10, iV e ~ 60, we find t r ~ 
10~ 25 seconds. In this case the maximum temperature of the Universe is 
~ 10 10 GeV. This is not very different to what one finds in simple inflationary 
models. 




(35) 




(36) 




(37) 
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As 4> rolls down the hill, one can check that leaves the slow-roll regime 
when e~ c< ^ exceeds 3/c 2 . At this point the constant term Vo in the potential 
becomes irrelevant and one can use the scaling solution for 0, all the way to 
the potential minimum. This is also the point at which density perturbations 
start to be generated via the ekpyrotic mechanism, while the Einstein frame 
scale factor a is still expanding. The Universe continues to expand slowly, 
but with a slowly decreasing Hubble constant, and finally enters contraction 
when when the density in the scalar field reaches zero, at a negative value of 
the potential energy. The ensuing contracting phase is accurately described 
by the scaling solution (|22|), in which a ~ (— t) p and = \fZp~jt, with t < 0, 
and t = being the time of the next big crunch. From the formulae (0) 
one finds di = (pai — ^Jp/3a )/t, which is greater than zero for p < |, since 
ao is greater than a±. Thus even when a is undergoing slow contraction, in 
the scaling era, the effect of the motion of is enough to make a\ expand 
throughout this phase. Matter residing on this brane would see continuous 
expansion all the way to the big crunch. The same argument shows that a-^ 
actually undergoes a small amount of contraction in the very much shorter 
scaling epoch of the expanding phase. 

6 Generation of Density Perturbations 

In the cyclic scenario, the period of exponential expansion occurring late in 
each cycle plays a key role in diluting the densities of matter, radiation and 
black holes to negligible levels, suppressing long wavelength perturbations 
and establishing a 'clean slate', namely a flat vacuous Universe in which all 
fields are in their quantum mechanical ground state. As the scalar field rolls 
down the potential in Eq. ([!]), entering the scaling solution in Eq. (^2]), the 
ekpyrotic mechanism for the generation of fluctuations derived in Refs. 1 
and [12] sets in and a scale invariant spectrum of adiabatic perturbations 
is thereby developed. Quantum fluctuations of the usual inflationary sort 
are also developed in the slow-roll quintessence epoch, but these are: (a) 
negligible in amplitude because Vo is tiny; and, (b) only excited on scales of 
order to and above in the contracting phase. These scales are shrunk only as 
(— t) p in the contracting, scaling solution, but then expanded as ta, t? and 

2 

tz in the kinetic dominated, radiation and matter eras in the big bang phase, 
which also lasts for a time of order to- Therefore, the modes amplified during 
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inflation are exponentially larger in wavelength than the Hubble radius scale 
in the next cycle by time of quintessence domination, which is the present 
epoch. 

Let us concentrate on the fluctuations produced via the ekpyrotic mecha- 
nism.0 Expanding the inhomogeneous fluctuations in the scalar field 5(p(t, x) = 
HzSfe^e x , we remove the damping term by setting 5<p^ = a~ 1 Xj:, to obtain 

4 = "* 2 XS + (j - V M cAj X % = ~{k 2 - 4) XS, (38) 

where primes denote conformal time derivatives and we have defined kp, 
the comoving 'freeze-out' wavenumber. Modes with k > kp oscillate with 
fixed amplitude, whereas those with k < kp are amplified. In the regime 
of interest kp grows monotonically so that shorter and shorter wavelengths 
progressively freeze out as the big crunch is approached. The physical scale 
at which modes freeze out is given by 

(39) 

As usual we adopt units where 8nG = 1, and denote proper time derivative 
with dot. In the era of quintessence domination when V dominates over 
Vj,,/,, the freeze-out scale is nearly constant, and comoving wavelengths 
are exponentially stretched beyond it. As begins to dominate however, 
Hubble damping becomes irrelevant, and the system approaches the scaling 
solution given in Eq. in which ~ — 2/t 2 , where t is the proper 

time to the big crunch. The freeze-out scale drops linearly with time to 
zero, as the scale factor is falling, like (— t)s. Therefore progressively shorter 
and shorter wavelength modes are frozen out and amplified, with waves of 
physical wavelength tp being frozen out at a time tp. 

An exponentially large band of comoving wavelengths is amplified and 
frozen in as (f) rolls from = down towards <p m in- Modes with all physical 
wavelengths from the microphysical scale t min , which could be not much 
larger than the Planck length, to the macroscopic scale to/c which is of order 
a tenth the present Hubble radius, aquire scale invariant perturbations. Once 
the perturbations are generated, their wavelength scales as (t min /tF) p in the 
collapsing phase. Then as crosses the potential well and races off to minus 
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infinity, the Einstein frame physical wavelength goes to zero. But this is not 
the relevant quantity to track, since we match the variables ao and a\ and 
therefore should match the physical wavelengths as measured by these scale 
factors. In the kinetic dominated phase, ao and a% are nearly constant, so in 
effect the physical wavelength of the modes are matched when crosses mm , 
in the contracting and expanding phases. Furthermore, the contracting and 
expanding solutions are nearly time-reverses of one another, until the time 
tDep computed above when the expanding solution deviates from scaling. 
Therefore one is effectively matching at tpep, from which one sees that the 
time tF at which perturbations on the current Hubble radius scale to were 
generated, is given by 

(40) 

where the bracketed factors are: (a) the contraction of the scale factor in the 
scaling solution, between the time £f at which the perturbations were gener- 
ated and the time toep at which the expanding solution departs from scaling; 
(b) the scaling back of the present comoving Hubble radius scale to the 
time of matter-domination £ m ; (c) the scaling back to the time of radiation- 
domination t r ; and (d) the scaling back to the time toe P using H 5 ~ constant, 
corresponding to kinetic domination in the expanding solution. 

From Eq. (|3(]), it follows that perturbations on the scale of the present 
Hubble radius were generated at a field value 



(41) 



Comparing with Eq. ([36|) for the resting value of the field 0c 1 , and the ex- 
pression 

2 

0mm ~ --In {ct /t min ) , (42) 

for the field value at the potential minimum, which follows from Eq. (|28|), 
one finds that 

0Gen - 0mm ~ ~^0mm + ^Tv^C ~ (fimin) , (43) 
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where the first term dominates. In other words, the fluctuations we see today 
were generated at a field value approximately half way between the zero and 
the minimum of V(<ft). 

7 Cyclic Solutions and Cyclic Attractors 

We have shown that a cyclic Universe solution exists provided we are allowed 
to pass through the Einstein-frame singularity according to the matching con- 
ditions elaborated in Section [|, Eqs. ( p^) and (pop. Specifically, we assumed 
that H 5 (out) = — (1 + x)Hr,(in) where \ is a non-negative constant, cor- 
responding to branes whose relative speed after collision is greater than or 
equal to the relative speed before collision. Our argument showed that, for 
each x > 0, there is a unique value of H 5 (out) that is perfectly cyclic. In the 
Appendix, we show that an increase in velocity is perfectly compatible with 
energy and momentum conservation in a collision between a positive and 
negative tension brane, provided a greater density of radiation is generated 
on the negative tension brane. (A similar outcome can occur through the 
coupling of (f> to the matter density, as discussed below Eq. fl25|), but we will 
only discuss the first effect for the purpose of simplicity.) 

In this Section, we wish to show that, under reasonable assumptions, the 
cyclic solution is a stable attractor, typically with a large basin of attraction. 
Without the attractor property, the cyclic model would seem fine-tuned and 
unstable. One could imagine that there would still be brane collisions and 
periods of contraction and expansion, but there would be no regularity or 
long-term predictability to the trajectories. If this were the case, fundamental 
physics would lose its power to explain the masses and couplings of elemen- 
tary particles. The masses and couplings depend on and other moduli 
fields. If there were no attractor solution, the precise trajectory of through 
cosmic history would depend on initial conditions and could not be derived 
from fundamental physics alone. In our proposal, the nature of the attrac- 
tor solution depends on microphysics at the bounce which is computable, in 
principle, from fundamental theory. Hence, masses and couplings of particles 
change during the course of cycle, but fundamental theory retains predictive 
power in determining the way they change and, specifically, their values at 
the current epoch. 

The essential feature for attractor behavior is the extended period of 
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accelerated expansion that damps the motion of 0. Let us consider how 
this works. Assuming x is fixed by microphysics, there is a value H^iput) 
which corresponds to the cyclic solution. Now, suppose the value of H 5 (out) 
exceeds H 5 (out). This means that the outgoing velocity exceeds the cyclic 
value and runs out farther on the plateau than in the cyclic case. Once the 
field stops, turns around, and quintessence-domination begins, the field is 
critically damped. By the time V(4>) falls to zero, the transient behavior of 
which depends on the initial value of H§ has damped away exponentially so 
that the field accurately tracks the slow-roll solution. Following the solution 
forwards, H§(in) at the next bounce is then exponentially close to what it 
would have been for the cyclic solution. By erasing memory of the initial 
conditions, the acceleration insures that H§(out) after the next bounce is 
very nearly H 5 (out). 

How many e-foldings of accelerated expansion are actually required to 
make the cyclic solution an attractor? If there is no epoch of accelerated 
expansion, perturbations will grow each cycle, becoming self-gravitating and 
non-linear so that no attractor will occur. A minimal requirement for ob- 
taining an attractor is that linear density perturbations grown during the 
matter era should be damped away during the subsequent exponential ex- 
pansion. This requires at least ln(10 5 )~ 10 e-foldings of exponential expan- 
sion. Equally, diluting the number density of baryons below one per Hubble 
volume is certainly over-kill in terms of ensuring an attractor, and this re- 
quires of order 60 e-foldings. In fact, as we discussed above, obtaining a far 
larger number of e-foldings is perfectly possible. 

To discuss the nature of the attractor solution, it is helpful to plot the 
trajectories of the system in the phase space given by the (H 5 , 0)-plane, 
shown in Figures |4] through |6|. Recall that H5 is proportional to 0; see 
Eq. ([0]). Figure |] illustrates the cyclic trajectory for the case where no 
radiation is generated at the bounce (x = 0) and the cycle is exactly time- 
symmetric. 

The phase space plot must always satisfy three properties. First, for a 
flat Universe, the Friedmann constraint equation H 2 = \p implies that the 
energy density p = (|0 2 + V(<f>)) must be positive. Without negative space 
curvature, the system is simply not allowed to explore negative energies. We 
show the classically excluded region as the shaded area in the Figures. In 
Figure (| where there is no radiation and V(4>) — > as — » 0, the excluded 
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region extends along H 5 = out to <fi — > — oo. In Figures || and |6] the 
shape of the excluded region is modified due to the presence of radiation. 
For example, the grey region pinches off on the left hand side for some finite 
value of 0. However, the effect is negligible for the trajectories considered in 
our discussion and so we show the same excluded region in the Figures as in 
the case of no radiation. 

The second property is that phase space trajectories are double-valued 
on the (i/5, 0)-plane. Given the scalar field and if 5, one may have either a 
contracting or an expanding Universe. We represent expanding trajectories 
as solid lines and the contracting trajectories as dashed lines. Two expanding 
trajectories are not allowed to cross, and neither are two contracting trajec- 
tories for the usual reasons that hold for particle trajectories on phase space. 
However, an expanding trajectory may certainly intersect a contracting tra- 
jectory. 

The final rule is that there are only two ways an expanding trajectory can 
turn into a contracting trajectory. If reversal occurs at finite <p h can only 
happen if the trajectory hits the forbidden zero density region (shaded), since 
p has to vanish if H is to pass smoothly through zero. The shaded region is 
analogous to the "egg" region described by Brustein and Veneziano." The 
second way in which contraction can turn into expansion is if the system runs 
off to = —00. Then, the 'bounce' described in Section [3] and Ref. || occurs. 

The trajectory shown in Figure f| is a cyclic solution (albeit not a very 
interesting one) in which no matter-radiation is produced at the bounce and 
the value of H 5 (out) is precisely equal to H 5 (in). The field rolls out in 
the expanding phase (solid line emanating from the upper left side of the 
Figure) to the value where V(<p) = and stops (the rightmost tip of the 
grey region). The total energy density is momentarily zero and expansion 
reverses to contraction. The field then rolls back to —00 (lower left side of 
the Figure). The expanding and contracting phases are exactly symmetrical. 

The time-scale for one cycle of this empty-Universe solution is easily esti- 
mated by noting that most of the time is spent near the zero of the potential, 
where it is rather shallow, and the scale factor a is nearly constant. There- 
fore, we can neglect gravity and calculate the period for one cycle in the case 
of the empty Universe: 
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Figure 4: The cyclic trajectory in the (H5, </>)-plane for the case where 
no matter and radiation are produced at the bounce (x = 0). The grey 
region which corresponds to negative energy density, is forbidden. The solid 
(dashed) line represents the trajectory during an expanding (contracting) 
phase. Expansion turns to contraction and vice versa when the trajectory 
hits the zero energy surface (the rightmost tip of grey region in this case). 
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for large c\(f) m i n \. For the parameters typical in the our examples, this corre- 
sponds to roughly one tenth of the current age of the Universe, or a billion 
years. 

In Figure [|, we consider the case where radiation is produced at the 
bounce, x > 0. If H 5 (out) is too low compared to the cyclic value, the 
trajectory encounters the zero density boundary (grey region) and reverses 
to contraction. Solid curves represent the expanding phase of the trajectory, 
and dashed lines represent the contracting phase of the trajectory. Let us 
call this an "undershoot" solution. H§ is only constant if <fi and a are both 
increasing, or both decreasing. However, if the Universe reverses when H§ is 
still positive, then the scalar field kinetic energy is blue shifted and and H§ 
is rapidly driven to more and more positive values. The trajectory flies off 
to large positive H 5 and (the upper boundary). 

As one increases H 5 (out), the behavior of the system changes. For suf- 
ficiently large H 5 , the system avoids the zero energy surface entirely during 
the period when is increasing (the bold solid and dashed trajectory in Fig- 
ure |5|). The field "overshoots" the negative region of the potential and lands 
on the positive plateau. Exponential expansion begins, followed by a very 
slow roll of back towards the potential zero, This period appears as a long, 
thin excursion on the right hand side of the figure. \H§\ is small because the 
field is rolling slowing in the quintessence-dominated phase). 

The cyclic attractor solution lies between these undershoot and overshoot 
trajectories. Figure ^ shows trajectories with initial values of H 5 (out) both 
above and below the cyclic value (the middle curve). Here we can study the 
stability of the cyclic solution. Let's first consider a trajectory with H 5 (out) 
larger than the value in the cyclic solution. This trajectory is indicated by 
(1) in the Figure. Clearly, it overshoots the cyclic trajectory and undergoes 
a longer period of exponential expansion (the long excursion to the right). 
During the slow-roll epoch, the difference between this trajectory and the 
cyclic one damps away until it is exponentially small. The trajectory en- 
counters the zero density surface very slightly later than the cyclic solution 
does, and, therefore, reverses and ends up with a very slightly smaller value 
of H^{out) than that in the cyclic trajectory. Similarly, one can see that 
starting the system in state (2) with a smaller value of H 5 (out) than that of 
the cyclic trajectory, the system will inflate less and reverse earlier, ending 
up with a larger value of H 5 (out) than that of the cyclic trajectory. This dis- 
cussion shows that the trajectory is stable and that memory of initial decays 
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Figure 5: Trajectories in the (H5, 0)-plane for the case where there is radia- 
tion. The solid (dashed) curves represent the trajectory during an expanding 
(contracting) phase. The thin lines illustrate undershoot solutions and the 
heavy line represents an overshoot solution. 
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exponentially after just one cycle. 

Here we implicitly assumed that x is a constant independent of H 5 (in), 
the incoming velocity. In the Appendix, we obtain an expression for x i n 
Eq. ([5(]) in terms of the matter-radiation energy densities created on the 
positive and negative tension branes. Assuming the energy density on the 
negative tension brane p_ is significantly greater than the energy density 
created on the positive tension brane, we have x oc p_. Assuming the col- 
lision between branes occurs at a low velocity so that one is not far from 
the adiabatic limit, p_ should decrease with decreasing H 5 m. (Note that the 
low-velocity assumption has been made throughout since it is required for the 
moduli approximation.) Hence, we would anticipate that x rises monotoni- 
cally as the incoming velocity increases. This effect can alter the trajectories 
and the precise basin of attraction in detail, but does not alter the conclusion 
that a large basin of attraction exists. This is assured by having a potential 
plateau or, more generally, a region of the potential in which slow-rolls 
with total energy comparable to the current dark energy density. 

Quantum effects are also unlikely to affect the attractor solution. We 
have shown that solution is stable under small perturbations and here the 
perturbations remain small since |Kmn| and Vq are small compared to the 
Planck scale. 

What of the trajectories in Figure |5], for example, that run away to large 
0? For these, it is important to understand what happens as grows more 
positive. One possibility is that the potential V(4>) diverges as — > oo. Our 
example for V((j)) has an infinite plateau, but, as discussed in Section |2], this 
is not a general requirement. If V(<p) grows sufficiently, will bounce back 
towards — oo. Alternatively, the same effect can occur if the theory includes 
massless fields that couple to the scale factor on the negative tension brane, 
a\. The Lagrangian density then includes a term a\ip 2 oc l/a\. Increasingly 
positive corresponds to shrinking a\. Hence, this dynamical term can also 
create a force that causes to bounce back. The net effect is that rattles 
back and forth along the potential, possibly following a chaotic orbit.011 
These effects could enhance the basin of attraction for the cyclic solution. 
That is, some of these trajectories which we ignored in our undershoot and 
overshoot treatment may eventually hit the plateau with low velocity, at 
which point they would become drawn to the attractor solution. 

Finally, let us emphasize that we have only considered the issue of sta- 
bility in the context of the very simplified model studied here, with a single 
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Figure 6: Trajectories in the (H^, 0)-plane showing the attractor nature of 
the cyclic solution (the middle trajectory). Path (1) is an overshoot solution 
that begins with slightly greater velocity (H 5 ) than the attractor, bounces 
off the grey zero-energy surface, and then has a contracting trajectory whose 
value of H 5 is smaller in magnitude. Path (2) is an undershoot solution which 
begins with slightly less velocity in the expanding phase than the cyclic value 
and ends in a contracting phase with H§ having a slightly greater magnitude. 
Following the next bounce, therefore, overshoot turns into undershoot and 
vice-versa. In either case, the deviation from the attractor value shrinks. 



40 



scalar field 0, and the matching conditions discussed in Section It would be 
very interesting to generalize this discussion to include other moduli, matter 
which couples in a nontrivial way to as discussed in Section |, and and also 
discrete degrees of freedom such as a quantized four-form field, which may 
change from cycle to cycle so that the system really explores moduli space. 
The existence or otherwise of an attractor could well be relevant to the deter- 
mination of the relative abundances of dark energy, dark matter, baryons and 
photons in the Universe, and also to the values of the fundamental constants 
of nature. 

8 Implications 

The strengths of the cyclic model are its simplicity, its efficient use of all 
of the dominant elements of the Universe and the fact that it is a com- 
plete description of all phases of cosmic evolution. This can be contrasted 
with inflationary cosmology, a highly appealing theoretical model in its own 
right. Inflationary cosmology focuses on a brief epoch when the Universe was 
10 -35 seconds old. The model relies on assumptions about how the Universe 
emerged from the cosmic singularity. One must postulate the existence of 
a phase of rapid cosmic acceleration at very high energies, for which their 
exists no direct proof. (In contrast, the cyclic model relies on low energy cos- 
mic acceleration that has been observed.) Subsequent cosmic events, such 
as the recent transition from matter-domination to dark energy domination 
and cosmic acceleration, appear to have no direct connection to inflationary 
theory. 

Because the cyclic model ties the past, present and future evolution of the 
Universe in a tight, cross-correlated way, it has surprising explanatory and 
predictive power. In the introduction to this paper, we noted a number of the 
most challenging questions of cosmology and fundamental physics. In this 
section, we consider each of these questions (and more) and briefly describe 
the insights the cyclic model provides concerning their answers. 
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8.1 Why is the Universe homogeneous, isotropic and 
flat? 

The Universe is made homogeneous and isotropic during the period of the 
preceding cycle when quintessence dominates and the Universe is undergoing 
slow cosmic acceleration. This ensures that the branes are flat and parallel 
as they begin to approach, collide, and emerge in a big bang. Inflation also 
relies on cosmic acceleration, but driven by very high vacuum energy which 
produces an acceleration that is nearly 10 100 times faster. 

8.2 How were density inhomogeneities generated? 

In the cyclic model, the observed inhomogeneities in the Universe are gener- 
ated during the contracting phase when the scale factor is nearly static and 
gravitational effects are weak. Consequently, as in the ekpyrotic scenario, a 
nearly scale invariant spectrum of adiabatic, gaussian energy density fluctu- 
ations is generated. However, because the expansion rate is negligible and 
gravitational effects are weak, the tensor (metric fluctuation) spectrum is 
blue with an exponentially tiny amplitude at long wavelengths. 

Fluctuations are also created during the quintessence dominated phase, 
just as they are during inflation. However, because the energy density during 
the accelerating phase is 100 orders of magnitude smaller than in inflation, the 
resulting fluctuation amplitude is exponentially smaller in the cyclic model. 
These fluctuations also have wavelengths that exceed the current Hubble 
horizon. Hence, they are observationally irrelevant. 

8.3 What is the role of dark energy and the current 
cosmic acceleration? 

Clearly, dark energy and the current cosmic acceleration play an essential 
role in the cyclic model both by reducing the entropy and black hole density 
of the previous cycle, and triggering the turnaround from an expanding to 
a contracting phase. (In all other cosmologies to date, including inflation, 
dark energy has no essential role.) 



42 



8.4 How old is the Universe? 



A truly cyclic Universe is clearly infinitely old in terms of cosmic time. As 
we have noted, the exact cyclic solution can also be an attractor. Hence, 
the cycling is stable. Consequently, one becomes insensitive to the initial 
conditions for the Universe as long as they were within the basin of attraction 
of the cyclic solution. We believe that within this framework, the problem of 
the initial conditions for the universe is significantly altered: as long as the 
universe has some nonzero probability for entering the cyclic solution, large 
regions of the universe maintain cyclic evolution for arbitrarily long periods 
of time. 

There is a possible objection to this argument, due to the fact that the four 
dimensional nonsingular brane spacetimes in our scenario are past geodesi- 
cally incomplete. As we have explained, for most of cosmic time they are 
well approximated by de Sitter spacetime, with a cosmological constant (or 
vacuum energy) close to the currently observed value. This nearly de Sitter 
spacetime is foliated by slices of constant scalar field <ft, which are nearly ge- 
ometrically flat. Matter is repeatedly generated on the slices with = — oo, 
in the rest frame defined by those slices. 

As one follows cosmic time t backwards, one must pass an infinite num- 
ber of these big crunch-big bang surfaces. However, even though the cosmic 
time tends to — oo, the proper time as measured along timelike geodesies 
running into the past generically is finite even as t tends to — oo. This may 
be seen as follows. Consider a particle with momentum P in the flat slicing. 
Its momentum blueshifts as a~ l as you follow it back in time. The geodesic 
becomes nearly null and the proper time measured along the geodesic con- 
verges even though t tends to — oo (this is the crux of the recent argument of 
Borde, Guth and Vilenkin that inflation is past geodesically incompleted) . 

In our scenario, however, all physical particles are created with finite mo- 
mentum in the flat slicing defined by (p. If we follow a particle present in 
today's universe back in time, most likely it was created on the last <ft — ~ oo 
surface. With an exponentially smaller probability, it could have been cre- 
ated on the penultimate <fi = — oo surface, and so on into the past. The 
probability that any observed particle originated on the t = — oo flat sur- 
face which is the boundary of the flat slicing of de Sitter spacetime, is zero. 
Therefore we do not attribute any physical significance to the past geodesic 
incompleteness of the spacetime metric in our scenario. In contrast, particle 
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production in standard inflationary models occurs on open slices a finite time 
ago. 

Even if there are no particles present which 'saw' the past boundary of 
the cyclic universe, one might object that the scenario implicitly requires 
a boundary condition in the infinite past. We do not think this is a strong 
objection. If the cyclic solution were begun within a finite region (for example 
a torus) of three dimensional space, it would grow exponentially with each 
cycle to an arbitrarily large size. After an arbitrarily long time, to any real 
observer the universe would appear to be infinite both in spatial extent, and 
in lifetime to the past. 

So, whilst the cyclic model still requires an initial condition, provided that 
state is within the basin of attraction of the cyclic solution, we are completely 
insensitive to its details. Any features of the initial state (the total size of 
the Universe, or any fluctuations about flatness or homogeneity), become 
exponentially diluted in each cycle and since the cyclic solution can repeat 
forever, are ultimately completely irrelevant to any observation. 

8.5 What is the ultimate fate of the Universe? 

The cycles can be continued to the infinite future, as well as the infinite past. 
Hence, the Universe endures forever. 

8.6 How big is the Universe? 

From the effective 4d point-of-view, the Universe oscillates between periods 
of expansion to periods of contraction down to a big crunch. However, from 
the brane world point-of-view, the Universe is always infinite in the sense that 
the branes always have infinite extent. The fact that the branes are spatially 
infinite means that it is possibly for the total entropy in the Universe to 
increase from cycle to cycle, and, at the same time, have the entropy density 
(in particular, the total entropy per Hubble horizon) become nearly zero prior 
to each bounce. 

8.7 What occurs at the big bang singularity? 

The cyclic model utilizes the ekpyrotic notion that the singularity corre- 
sponds to the collision and bounce of two outer orbifold branes in a manner 
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that is continuous and well-behaved. The singularity is not a place where 
energy and curvature diverge and time begins. Rather, formulated in ap- 
propriate fields and coordinates, the singularity is a smooth, finite transition 
from a contracting phase heading towards a big crunch and a big bang evolv- 
ing into an expanding Universe. 

8.8 What determines the arrow of time? 

Since the Universe is cyclic, it may appear that there is no well-defined means 
of determining the arrow of time. Indeed, for a local observer, there is no 
clear means of doing so. 

From the global perspective, though, there is a clear means of determining 
forward from backward in time. First, one of the boundary branes is forever 
expanding in the "forward" time direction in the cyclic model. The other 
brane is expanding except for brief intervals of contraction, but, averaged 
over a cycle, the net effect is expansion. The rate changes from phase to 
phase, as well as the separation. In the contraction phase, the branes them- 
selves stretch at a rate that is slow but their separation rapidly decreases. In 
the radiation, matter, and quintessence dominated phases, the branes stretch 
significantly, but their separation remains fixed. During this period, the en- 
tropy created during the previous cycle is spread out exponentially, reducing 
the degrees of freedom per horizon to nearly zero. 

8.9 Why is the cosmological constant so small? 

The cyclic model provides a fascinating new outlook on this vexing problem. 
Historically, the problem is assumed to mean that one must explain why the 
vacuum energy of the ground state is zero. 

In the cyclic model, the vacuum energy of the ground state is not zero. 
It is negative and its magnitude is large, as is obvious from Figure 1. If 
the Universe begins in the ground state, the negative cosmological constant 
will cause rapid recollapse, as expected for an anti-de Sitter phase. In the 
cyclic scenario, though, we have shown how to arrange conditions where the 
Universe avoids the ground state. Instead, the Universe hovers from cycle to 
cycle above the ground state bouncing from one side of the potential well to 
the other but spending most time on the positive energy side. The branes 
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are moving too rapidly whenever the separation corresponds to the potential 
minimum. 

There remains the important challenge of explaining why the the cur- 
rent potential energy is so small. The value depends on both the shape of 
the potential curve and the precise transfer of energy and momentum at the 
bounce. Perhaps explaining the value will be an issue as knotty as the cos- 
mological constant problem, or perhaps the conditions will prove easier to 
satisfy. What is certain, though, is that the problem is shifted from tuning a 
vacuum energy, and this provides an opportunity for new kinds of solutions. 



8.10 Equation-of-state of dark energy 

The equation-of-state of the dark energy, w, is the ratio of the pressure to 
the energy density of 0, (|0 2 — V)/{\4> 2 + V). In Section |, we discussed the 
evolution of <ft i n the radiation, matter and quintessence dominated epochs. 
The generic result is that evolution of <fi in the positive direction halts and 
the field begins to roll back towards — oo in the radiation-dominated epoch. 
At the turn-around, w = —1 since the kinetic energy is zero. As the field rolls 
back and its kinetic energy increases, w increases. Hence, the generic result 
is that w is close to -1 today and increasing. Conceivably, cosmological 
observations could detect this prediction. Tracker models of quintessence, 
some of the best-motivated alternatives, have the opposite trend: w is near 
—0.8 or so today and decreasing towards —1.0 Other models, such as k- 
essence, have the same trend as found in the cyclic model.il 



8.11 Implications for Supersymmetry and Superstrings 

The cyclic model imposes different constraints on fundamental physics com- 
pared to previous cosmological models. As an example, consider the problem 
of designing supergravity potentials. The potentials are constructed from a 
superpotential W according to the prescription: 
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state is super symmetric, D{W = 0, the first term is zero. In general, unless 
W is zero for precisely the same values for which DiW = 0, the minimum 
has a negative cosmological constant. In the past, this type of model would 
have been ruled unacceptable. The possibility of a cosmology in which the 
Universe hovers over the ground state in a state of zero or positive energy 
revives these models and alters constraints on model building. 

An obvious but important implication is that supersymmetry breaking 
can be achieved without having spontaneous symmetry breaking in the ground 
state. In this scenario, it suffices if the Universe hovers in the radiation, mat- 
ter and quintessence dominated epochs at some state far above the ground 
state in energy and that the supersymmetry is broken by the appropriate 
amount in the hovering state, where the radiation, matter and quintessence 
dominated phases occur. These considerations have a significant impact on 
the design of phenomenological supersymmetric models. 

One other requirement/prediction of the cyclic scenario (and the ekpy- 
rotic models in general) is that the branes move in a space-time with co- 
dimension one. The constraint derives from having a bounce that produces 
a smooth transition from contraction to expansion. As argued by Khoury 
et ai, the geometry is flat arbitrarily close to the bounce provided there is 
one extra dimension only. Hence, brane world scenarios based on theories 
like that of Horava and Witten are acceptable, but large extra-dimensional 
models relying on having co-dimension two or greater are problematic. 

8.12 Hoyle's Revenge? 

Within each cycle, there is a sequence of kinetic energy, radiation, matter and 
quintessence dominated phases of evolution that are in accord with the stan- 
dard big bang cosmology. However, averaged over many cycles, the model can 
be viewed as a remarkable re-incarnation of Fred Hoyle's steady state model 
of the Universe. Most of the cycle is spent in a phase with nearly constant 
energy density, as in the steady state picture. Indeed Hoyle's C-field that 
was introduced to provide a constant supply of matter (and a preferred rest 
frame) is replaced by our scalar field <f>, which defines a preferred time slicing 
and generates matter repeatedly at each bounce, restoring the Universe to a 
state of high temperature and matter density. In Hoyle's steady state model, 
every flat spatial slice was statistically identical. Here the slices are iden- 
tical only when separated by one period, so we have a discrete rather than 
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continuous time translation symmetry. Nevertheless when coarse grained 
over large time spans, the structure is similar to that proposed in the steady 
state Universe. Global properties of the cyclic cosmology will be discussed 
elsewhere.il 
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9 Appendix: Matching H$ across the Bounce 

In this appendix we discuss the matching condition needed to determine 
Hz,{out) in terms of H§{in). We shall assume that all other extra dimensions 
and moduli are fixed, and the bulk space-time between the branes settles 
down to a static state after the collision. (In the simplest brane world models, 
there is a Birkhoff theorem which ensures that there is a coordinate system 
in which the bulk metric is static in between the branes). We shall take the 
densities of radiation on the branes after collision as being given. By imposing 
Israel matching in both initial and final states, as well as conservation of total 
energy and momentum, we shall be able to completely fix the state of the 
outgoing branes and in particular the expansion rate of the extra dimension 
H^out), in terms of if 5 (in). A more complete discussion of this method will 
be presented in Ref. 

The idea is to treat the brane collision as a short- distance phenomenon. 
The warp factor may be treated as linear between the branes as they approach 
or recede. Linearity plus Z2 symmetry ensures that the kinks in the warp 
factors are equal in magnitude and opposite in sign. Israel matching relates 
the kink magnitudes to the densities and speeds of the branes, yielding the 
relations we use below. 

The initial state of empty branes with tensions T and — T, and with 
corresponding velocities v + < and > (measured in the frame in which 
the bulk is static) obeys 
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The first equation follows from Israel matching on the two branes as the 
approach, and equating the kinks in the brane scale factors. The second and 
third equations are the definitions of the total energy and momentum. The 
three equations (pE6|) imply that the incoming, empty state has v + = — 1>_, 
E to t = and that the total momentum is 

°„ = , TL "^ < 0, (47) 



{\ - \(LH 5 (m))- 



where we identify v + — V- with the contraction speed of the fifth dimension, 
\LH 5 (in)\. 

The corresponding equations for the outgoing state are easily obtained, 
by replacing T with T + p + = T + for the positive tension brane, and — T 
with — T + p_ = — T_ for the negative tension brane, assuming the densities 
of radiation produced at the collision on each brane, p + and p_ respectively, 
are given from a microphysical calculation, and are both positive. 

We now wish to apply energy and momentum conservation, and Israel 
matching to the final state. The only subtlety is that the (t, y) frame in 
which the bulk is static is not necessarily the same frame in the final state 
as it was in the initial state, so one should boost the initial two-momentum 
{E t ot, Ptot) with a velocity V and then apply the Israel constraints and energy- 
momentum conservation equations in the new boosted frame. The latter 
provide three equations for the three unknowns in the final state, namely 
v + (out), v^(out) and V. Writing v±(out) = tanh(6*±), where 9± are the 
associated rapidities, one obtains two solutions 

sinh0 + = -^-{\P tot \ + \P tot \-\T 2 + -T 2 



2T_ 

smh#_ = ±^-{\P tot \-\ Ptot \-\Tl-T 2 )), (48) 

where T + = T + p + , T_ = T — p_ with p + and p_ the densities of radiation on 
the positive and negative tension branes respectively, after collision. Both p + 
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and p- are assumed to be positive. In the first solution, with signs ( — h), the 
velocities of the positive and negative tension branes are the same after the 

collision as they were before it. In the second, with signs ( ), the positive 

tension brane continues in the negative y direction but the negative tension 
brane is also moving in the negative y direction. 
The corresponding values for v±(out) and V are 



v + (out) 
V- (out) 

V 

or 



Ptat\ + \Pt 



tot 



T 2 



pL+2(n 



T 2 )+P t -J(Tl-T^ 



\P 



tot 



\p 



tot I 



1-1 



(T* 



T 2 



pL+2(n 



T 2 



PfJ(T 2 -T 2 ) 2 



P? 



tot 



2 (T 2 +T 2 ) + P to 2 (T 2 -T 2 ) 2 



\Ptot\{T 2 + + T 2 )/(T 2 - T 2 ) + \Ptot\-\Tl - T 2 ) ' 



pL+nn 



T 2 ) + p t -j(Ti-T 2 y 



iPotl + lPtotl-HTl + Tl 



(49) 



where the first solution for V holds for the ( — h) case, and the second for the 
( — ) case. 

We are interested in the relative speed of the branes in the outgoing 
state, since that gives the expansion rate of the extra dimension, —v + (out) + 
V-(out) = LH 5 (out), compared to their relative speed — 2v + = —LH 5 (in) in 
the incoming state. We find in the ( — h) solution, 



H 5 (out) 



H 5 (in) 
and in the 



v + (out) — v_(out) 
2vl 



solution 



p2 

r tot 



AT 2 



\p t 2 ot + 2(T 2 +T 2 )+P t -J(T 2 -T_ 



2\2 : 



(50) 



H 5 (out) 



HMn) 



(T 2 - T 2 

r tot 



Pit + 4T 2 

\ P 2 ot + 2(T 2 + T 2 ) + P t -J(T 2 - T 2 ) 2 • 



(51) 



with P to t given by (^) in both cases. 

At this point we need to consider how the densities of radiation p + and p_ 
depend on the relative speed of approach of the branes. At very low speeds, 
\LH 5 (in)\ « 1, one expects the outer brane collision to be nearly adiabatic 
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and an exponentially small amount of radiation to be produced. The ( — h) 
solution has the speeds of both branes nearly equal before and after collision: 

we assume that it is this solution, rather than the ( ) solution which is 

realised in this low velocity limit. 

As \LH 5 (in)\ is increased, we expect p + and p_ to grow. Now, if we 
consider p + and p_ to be both << P tot « T, then the second term in 
the denominator dominates. If more radiation is produced on the negative 
tension brane, p_ > p + , then \H 5 (out) / H 5 (in)\ = (1 + x) ~ (1 + (p_ — 
p + )/2T) and so x is small and positive. This is the condition noted in the 
text, necessary to obtain cyclic behavior. Conceivably, the brane tension 
can change from T to T" = T — t at collision. Then, we obtain (1 + x) ~ 
(l + (p_-p+ + 2t)/2T) 

For the ( — h) solution, we can straightforwardly determine an upper limit 
for \H 5 (out) I H 5 (in)\ = (1 + x)- Consider, for example, the case there the 
brane tension is unchanged at collision, t = 0. The expression in ([5(]) gives 
\H 5 (out) / H 5 (in)\ as a function of T+, TL and P to t- It is greatest, at fixed T_ 
and Ptot, when T + = T, its smallest value. For P 2 ot < T 2 , it is maximized 
for T 2 =T 2 - P 2 ot , and equal to y/l + P t 2 ot /{4T 2 ) when equality holds. For 
Ptot > T 2 , it is maximized when T_ = 0, its smallest value, and P 2 ot = 2T 2 , 
when it is equal to y|. This is more than enough for us to obtain the small 
values of x needed to make the cyclic scenario work. A reduction in brane 
tension at collions t > further increases the maximal value of the ratio. 
To obtain cyclic behavior, we need x to be constant from bounce to bounce. 
That is, compared to the tension before collision, the fractional change in 
tension and the fractional production of radiation must be constant. 

We shall not consider the ( ) solution in detail, except to note that in 

the small P to t limit it allows an arbitrarily large value for \H 5 (out)/ H 5 (in)\, 
which seems unphysical. 

Let us re-iterate that, there are many caveats attached to this calculation. 
We have not calculated p + and p_ and have left these as parameters. We 
have set the brane tensions T to be equal before and after the collision, 
and have neglected possible bulk excitations, treating the bulk spacetime as 
static both before and after collision. We have ignored matter couplings to 
bulk scalars, and ignored the possible dynamical evolution of additional extra 
dimensions. Nevertheless we think it encouraging that the unusual behavior 
of matter bound to a negative tension brane allows H^(out) / H^,(in) to be 
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slightly greater than unity, which is what we need for cyclic behavior. 
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